Antiphase dynamics in a multimode semiconductor laser with optical injection 
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A detailed experimental study of antiphase dynamics in a two-mode semiconductor laser with 
optical injection is presented. The device is a specially designed Fabry-Perot laser that supports 
two primary modes with a THz frequency spacing. Injection in one of the primary modes of the 
device leads to a rich variety of single and two-mode dynamical scenarios, which are reproduced with 
remarkable accuracy by a four dimensional rate equation model. Numerical bifurcation analysis 
reveals the importance of torus bifurcations in mediating transitions to antiphase dynamics and of 
saddle-node of limit cycle bifurcations in switching of the dynamics between single and two-mode 
regimes. 

PACS numbers: 42.55.Px, 42.65.Sf, 05.45.-a 



I. INTRODUCTION 

The semiconductor laser with optical injection is a con- 
ceptually simple system that nevertheless has complex 
nonlinear dynamical properties. Phenomena of interest 
that have been observed experimentally include subhar- 
monic resonance [l[ , period doubling route to chaos [1, Q 
and bistability of wave-mixing and injection locked states 
[3, @ ■ Use of numerical continuation tools has enabled 
diverse dynamical states of the single mode edge-emitting 
system with optical injection to be linked together in a 
consistent way ;6;] . More recently, a global picture of the 
complex dynamics of polarization modes in VCSELs with 
orthogonal optical injection has also begun to emerge 

For both fundamental and applied reasons, interest 
in multimode dynamics of semiconductor lasers has in- 
creased. This follows from their increased dynamical 
complexity, as well as possible applications in optical sig- 
nal processing [1, Q. Multiple oscillating modes allow 
for antiphase behaviour, where mutual coupling of indi- 
vidual modal degrees of freedom leads to anticorrelated 
intensity dynamics 0, [H [ll [S [3 ■ 

Recently, we have developed a specially engineered 
two-mode diode laser with a large (THz) primary mode 
spacing [l^ . The device can be biased such that it oscil- 
lates on two modes simultaneously with the same average 
power level. Unlike the case of polarization modes in VC- 
SELs, the mode spacing is in the highly non-degenerate 
regime, allowing us to develop studies of the multimode 
nonlinear dynamics in the minimal possible two mode 
system. These studies provide a means to test models of 
the multimode diode laser and to reveal new bifurcation 
structure that occurs beyond the single mode approxi- 
mation jlGj . 

A number of approaches have been developed to mod- 
eling the multimode dynamics of edge-emitting semicon- 
ductor lasers. Partial differential equations will generally 
reproduce experimental results, but such an approach 
suffers from computational difficulty and the reduced 
level of physical insight provided [l3, Ell. Rate equa- 
tion approaches have proven very successful in the case 



of modeling the dynamical response of single mode semi- 
conductor lasers. In the case of multimode dynamics, 
however, a number of different models have been devel- 
oped that differ in how the d yna mics of the carrier den- 
sity are described [H, [H, [13, IHI . 

Here we provide an experimental and theoretical study 
of antiphase dynamics when one of the two primary 
modes of the device is optically injected. Among the 
rich variety of multimode dynamical phenomena we have 
found limit-cycle, quasiperiodic and chaotic antiphase 
dynamics as well as regions where the optical power is 
largely switched to the uninjected mode. As we will show, 
all of these dynamics, as well as the familiar injection 
locking on a single mode, can be found by performing a 
sweep over a large detuning range of the injected field 
for a certain fixed injected field intensity. As the cen- 
tral theme of this study, we provide a detailed descrip- 
tion of regions of quasiperiodic antiphase dynamics using 
spectrally resolved power spectral measurements and in- 
tensity time traces. We then demonstrate that a four 
dimensional rate equation model reproduces the overall 
experimental picture very well. In particular, we show 
how regions of quasiperiodic antiphase dynamics that we 
have observed experimentally are reproduced in numeri- 
cal simulations and that new bifurcation structure, which 
is particular to the two mode system considered, governs 
the appearance of these dynamics. 

This paper is organised as follows. In section|TTl we de- 
scribe our device and experimental set-up and we provide 
optical and power spectral data that illustrate the variety 
of dynamical scenarios that can be observed as a function 
of frequency detuning at a fixed injected field strength. 
We also provide intensity time traces and power spec- 
tral data characterizing two regions of quasiperiodic an- 
tiphase dynamics. In section [TTTl we introduce our model 
equations that describe the two-mode injected system, 
showing how the overall experimental picture as well as 
specific examples of antiphase dynamics are qualitatively 
reproduced. Finally, in section [TVl we describe the bifur- 
cation structure that leads to the antiphase dynamics 
observed. 
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FIG. 1: (Color online) Optical spectra of the two-colour de- 
vice as the frequency detuning, /S.lu is varied at a fixed in- 
jection. Left: uninjected mode, u\. Right; Injected mode, 
fi- 



ll. EXPERIMENT 

For the primary mode spacing of four fundamental FP 
modes (480 GHz), we observe simultaneous lasing of the 
two primary modes. We adjust the device current in or- 
der that the time averaged optical power in each primary 
mode of the free running laser is equal. Details of the de- 
sign and free running lasing characteristics of the device 
we consider can be found in Ref p^ . 

For our experiment, we use the two-colour laser in a 
master-slave configuration, where the master laser is a 
tunable laser with < 100 kHz linewidth. Choosing a 
fixed injected power level, we vary the detuning, Atj, 
of the injected field from the long wavelength mode of 
the device, V2- The optical spectra of the two primary 
modes are shown as a function of the frequency detun- 
ing of the injected signal in Fig. [T] for a detuning range 
of -9 GHz < Aw < 14 GHz. Despite the fact that we 
vary only the frequency detuning. Fig. [T] indicates that 
we are dealing with a rich dynamical scenario; regions of 
multiwavelength dynamics are interspersed with regions 
of single or nearly-single mode dynamics, where the op- 
tical power can be concentrated in either of the primary 
modes. In sections IIIII and IIVI we will show that a min- 
imal extension of the single mode rate equation model 
by a single dynamical variable is sufficient to account for 
the presence of the uninjected mode and reproduces the 
experimental data with remarkable accuracy. 

For this injected power level the power spectral density 
of the total intensity is plotted in Fig. [21 By considering 
the spectrum of the total intensity we naturally exclude 
the antiphase frequency components that are present in 
each region of two mode dynamics. From Figs [1] and [21 
the regions of stable and unstable injection locking where 
the uninjected mode is off can be identified near zero de- 
tuning. In the injection locking region, a sharp single 
peak at the frequency of the injected field is present in 
the optical spectrum at long wavelength. In the region of 
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FIG. 2: (Color online) Power spectral density of the total 
intensity as a function of frequency detuning. The injected 
power level is the same as in Fig. [T] 

unstable locking, which extends from approximately zero 
detuning to -2.5 GHz, a single frequency that is associ- 
ated with undamped relaxation oscillations is visible in 
Fig. [21 This frequency is almost constant at c. 5.5 GHz, 
which we identify as the relaxation oscillation frequency 
of the device. 

At large positive frequency detuning of 14 GHz, we find 
a broad region of nearly single mode dynamics, where 
the optical power is concentrated in the injected mode. 
In this region the injected field and two weaker symmet- 
rically placed components can be clearly resolved in the 
optical spectrum at long wavelength. Examining Fig. [2l 
we find that in the corresponding region only the detun- 
ing frequency is present. This is then a region of wave 
mixing or beating between the injected field and the long 
wavelength primary mode. 

As the detuning is decreased, a period doubling bifur- 
cation occurs and we enter a region where a new fre- 
quency equal to one-half of the detuning frequency is 
present in Fig. [21 As indicated in Fig. [Tl within this 
broad region we find region (i) of two mode dynamics, 
where the optical spectra of both modes are are consid- 
erably broadened. A representative set of power spectral 
densities and time traces of the total intensity and the in- 
dividual modal intensities in region (i) are shown in Fig. 
[21 In this figure the detuning is 10 GHz. As expected the 
power spectrum of the total intensity shown in Fig [31 (b) 
comprises two broadened peaks which include distinctive 
shoulders. 

In the power spectra of the intensities of the individual 
modes shown in Fig. [31 (d) and (f) we can see a strong 
antiphase frequency component with a distinct peak at 
c. 400 MHz. We see also that the features associated 
with the injection frequency and its subharmonic in each 
of these figures are broadened further because of the for- 
mation of satellites due to mixing with the low frequency 
antiphase component. Note also that the detuning of the 
injected signal is the dominant frequency in the dynam- 
ics of the injected mode [Fig. [31 (f )] but the subharmonic 
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FIG. 3: Experimentally measured antiphase dynamics, region 
(i) . The frequency detuning is 10 GHz and the injected power 
level is the same as in Fig. [1] Left panels: Intensity time 
traces. Right panels: Power spectral densities, (a) and (b): 
Total Intensity, (c) and (d): uninjected mode, (e) and (f): 
Injected mode. 



FIG. 4: Experimentally measured antiphase dynamics, region 
(ii). The frequency detuning is 4 GHz and the injected power 
level is the same as in Fig. [T] Left panels: Intensity time 
traces. Right panels: Power spectral densities, (a) and (b): 
Total Intensity, (c) and (d): uninjected mode, (e) and (f): 
Injected mode. 



near the relaxation oscillation frequency dominates in the 
dynamics of the uninjected mode [Fig. [3] (d)]. 

We have confirmed that the low antiphase frequency is 
an independent frequency component that is not a linear 
combination of the other frequencies present. The low 
frequency is clearly visible as an envelope modulation of 
the intensity in the time traces of the individual modes 
shown in Fig. [3] (c) and (e) . These antiphase time traces 
have a distinctive sawtooth structure that results from 
the growing or decaying oscillations that form each pulse. 
They are remarkably similar to the so-called regular pulse 
packages that can be observed in semiconductor lasers 
with optical feedback in the short external cavity regime 
[2^ . However, in the case of optical injection, the fast 
time scale is determined by the relaxation oscillation and 
its harmonics, rather than the external cavity frequency. 

As the detuning is decreased further a reverse period- 
doubling occurs and a simple single mode wave mixing 
region is found. Adjacent to this region, we next find a 
second broad region of two mode dynamics, which has a 
lower boundary at zero detuning. For reasons that will 
become clear, we subdivide this region into two, and label 
the sub-regions as regions (ii) and (iii). Broadly speak- 
ing, region (ii) is found once the detuning is less than 
the relaxation oscillation frequency. Note that across the 
boundary between the single mode wave mixing region 
and region (ii), the power spectrum of the total intensity 
briefly becomes strongly broadened. Then, as the detun- 
ing is decreased further, the dynamics become simpler 
with the relaxation oscillation frequency, and its har- 
monics, largely determining the dynamics of the total 
intensity. 

A representative set of power spectral densities and 
time traces in region (ii) from inside this boundary are 
shown in Fig. 01 The dynamics of the total intensity 



is now largely determined by the relaxation oscillation 
and its harmonic although weak and broadened satellite 
features are also visible [cf. Fig. [4] (b)] . 

From Fig. [4] (d) and (f), we see that the dynamics 
of the individual modes are determined by a series of 
frequencies. In particular, there is a strong low frequency 
component present in the individual dynamics of both 
modes that is almost completely in antiphase. As was 
the case in region (i), this low frequency component is 
responsible for the appearance of a cascade of frequencies 
that are centered at the relaxation oscillation peak and 
its harmonic. 

As the detuning is decreased further, the frequency 
difference between the relaxation oscillation and the de- 
tuning approaches half of the relaxation oscillation fre- 
quency. In this region a weak period two region of the 
total intensity is found. For detunings less than approx- 
imately one-half the relaxation oscillation frequency we 
find coupled chaotic dynamics that are characterised by 
broad-band spectra. 

For negative detunings, the single mode locking region 
is found, which extends to a detuning of approximately 
-5 GHz. At the boundary of the stable locking region, in- 
side of which there is no structure in the power spectrum, 
we find region (iv), where the spectrum is again strongly 
broadened. Decreasing the detuning further, the spec- 
trum is almost structureless. Here, a weak feature at the 
detuning frequency can be seen. This corresponds to the 
transfer of the optical power to the uninjected mode. 



III. MODELING OF THE DEVICE RESPONSE 

We have adapted the well known model of a single 
mode laser [l^, [2^ with optical injection to account for 
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the presence of a second lasing mode. The system of 
equations in normahzed units may be written as follows 

^l = i(l + ^a)(.9l(2n+l)-l)Sl (1) 

E2 = [5(1 + «a)(52(2n + 1) - 1) - iAtu] E2 + K (2) 
rn = P-n-(l + 2n)^^5„l£;,„|2 (3) 

where the nonlinear modal gain is 

9m^9'^°^(l + eY,f^^n\E^\^ . (4) 

Here Ei and E2 are the normalized complex electric 
field amplitudes of the modes and n is the normalized 
excess carrier density. The bifurcation parameters are 
the normalized injected field strength K and the an- 
gular frequency detuning Aw. Further parameters are 
the phase-amplitude coupling a, the product of the car- 
rier lifetime and the cavity decay rate T, the normal- 
ized pump current P, and the linear modal gain . In 
our numerical simulation we used the values a = 2.6, 
= 0.00125, P = 0.5, (twice threshold), and g^Jn' = 1- 
Then the value of the relaxation oscillation frequency is 
tjRO = \J'^P jP ^ 5.5 GHz, in agreement with experi- 
ment. The cross and self saturation are determined by 
e/3mri and we use the values e = 0.01, /3i2 = /32i = 2/3 
and = /?22 = Ij which is consistent with the stability 
of the two mode solution in the free running laser. 

Note that although we have provided a complex equa- 
tion for the field E\ , the phase of E\ is in fact decoupled 
leading to a four-dimensional system of equations where 
only the intensity of the uninjected mode influences the 
dynamics. In addition, the single mode dynamics is con- 
tained within the invariant submanifold {E\ — 0) in these 
equations. Our model equations are therefore a minimal 
extension of the (three dimensional) single mode case and 
follow from the fact that the primary mode spacing in the 
device considered is in the highly non-degenerate regime. 

This system only considers a single averaged carrier 
density and a general cross and self saturation of the gain. 
More complex models of the multiwavelength dynamics 
are available that are derived by considering additional 
Fourier components of the carrier density. One example 
includes the effect of static spatial hole burning explicitly 
leading to a system with two carrier density equations, 
one defined for each field [2l|. However, we have found 
that these equations lead to some unphysical results in 
the injected system, such as a stable injection locked re- 
gion where the uninjected mode is not suppressed. 

Another modeling approach [ij, [2^ includes an asym- 
metric term that is associated with a dynamic grating 
formed in the carrier density profile. This term effectively 
represents the mutual injection or wave mixing interac- 
tion of the primary modes and has been shown to play 
a role in the fre e running switching dynamics of Fabry- 
Perot lasers [1J|. This grating is also responsible for the 
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FIG. 5: (Color online) Power spectral density of the total field 
intensity + |-B2l^ obtained from numerical integration of 
Eqs. Il}- ([3]) as a function of the detuning. The injected field 
strength is if = 0.008. 

formation of four-wave mixing sidebands in the free run- 
ning two-mode laser considered here [151] . However, we 
have found equivalent experimental results for injection 
at both the long and short wavelength primary modes, 
suggesting that this grating does not significantly affect 
the dynamics. We attribute this to the very large sep- 
aration of the primary modes of the device (480 GHz), 
which means that the dynamic grating is weakly devel- 
oped, with large sidebands observed because of enhance- 
ment by the Fabry-Perot cavity. 

Using the model equations H])-®, in Fig. [5] we plot 
the power spectral density of the total field intensity 
\E\? + |i^2p for fixed K and for detuning -9 GHz < 
Aw < 14 GHz. The value of K is 0.008 which was chosen 
as it provided the best possible agreement with experi- 
ment over the whole range of frequency detuning. When 
compared with the experimental result of Fig. [21 we find 
that the level of qualitative agreement is excellent. 

To further illustrate that the individual modal dynam- 
ics are also well reproduced, in Fig. [6|we plot the local 
extrema of the numerically calculated field intensities, 
\E\\^ and |i?2p, over the same detuning range. 

If we examine the structure of Fig. [BJ we can immedi- 
ately identify two regions where the uninjected mode is 
off and that the position and extent of these regions is in 
good qualitative agreement with the experimental results 
shown in Fig. [T] Thus not only do the model equations 
reproduce the dynamics of the total intensity, they also 
provide an accurate picture of the regimes of single and 
multimode dynamics of the system. Note that this agree- 
ment extends to the model reproducing the strong sup- 
pression of the injected mode at large negative detunings 
that is observed experimentally. 

In Fig. [HI vertical arrows indicate the regions (i) to 
(iv) of two mode dynamics that we have identified exper- 
imentally. For detailed comparison with experiment, we 
now provide the numerical intensity time traces and cor- 
responding power spectra for the two regions of antiphase 
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FIG. 6: Local extrema of the field intensities (a) \Ei\^ and 
(b) I-E2P obtained from numerical integration of Eqs. ((T|- 
((3]) as a function of the detuning. The injected field strength 
is .K" = 0.008. Upper panel: uninjected field. Lower panel: 
Injected field. 
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FIG. 7: Numerically computed antiphase dynamics, region 
(i). The frequency detuning is 10.0 GHz. The injected field 
strength is K = 0.008. Left panels: Intensity time traces. 
Right panels: Power spectral densities, (a) and (b): Total 
Intensity, (c) and (d): uninjected mode, (e) and (f): Injected 
mode. 



dynamics labelled (i) and (ii). 

In Fig. [3 we show these data for the total intensity 
and for the intensities of each of the two modes at a de- 
tuning of 10.0 GHz, which is located in region (i). What 
is most striking is the reproduction of the sawtooth pulse 
package structure that was seen experimentally. The en- 
velope period of the pulse packages is seen to be c. 300 
MHz, in broad agreement with experiment. The cascade 
of equally spaced frequencies that begins at this enve- 
lope frequency are also almost completely in antiphase. 
The power spectrum of the total intensity is also seen 
to comprise two broad cascades formed by the envelope 
frequency, one centered at the detuning frequency and 
another at the relaxation oscillation frequency. However, 
the number of these peaks is much larger than was found 
experimentally in the case of the individual modes and, 
in fact, the cascade structure that is seen in the numer- 
ical data for the total intensity is poorly resolved in the 
experiment. 
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FIG. 8: Numerically computed antiphase dynamics, region 
(ii). The frequency detuning is 3.4 GHz. The injected field 
strength is K = 0.008. Left panels: Intensity time traces. 
Right panels: Power spectral densities, (a) and (b): Total 
Intensity, (c) and (d): uninjected mode, (e) and (f): Injected 
mode. 



In Fig. [51 we show these data for the total intensity 
and for the intensities of each of the two modes at a de- 
tuning of 3.4 GHz, which is located in region (ii). We 
have had to chose a slightly smaller value for the detun- 
ing frequency for this simulation because of the discrep- 
ancy that exists between the boundary of the central two 
mode region that was found numerically and in the ex- 
periment. Nevertheless, the qualitative agreement with 
the experimental time traces and power spectra is again 
very good in this region. We note in particular the ex- 
perimental intensity waveforms of the individual modes 
that are peaked in the case of the uninjected mode and 
rounded for the injected mode are well reproduced. 



IV. BIFURCATION ANALYSIS 

The bifurcations evident from Fig. [5] and Fig. [6] can 
be understood from the global bifurcation diagram in 
the Alli vs. K plane as shown in Fig. [9l Let us focus 
on the vertical line at injection strength K = 0.008, 
which was used in Fig. [5] and Fig. [6] and identified to 
best fit the experimental data shown in Figs. [T] and [D At 
large Au = 15 GHz we are in a two mode wave mixing 
region, with the uninjected mode vi largely suppressed 
and weakly modulated. This agrees with the experimen- 
tal observation of a weak but non-vanishing uninjected 
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FIG. 9: (Color online) Bifurcation diagram in the Alu vs. 
K plane of Eq. ©-([S]) calculated using the numerical con- 
tinuation tool AUTO-07p 24]. Solid and dashed lines denote 
bifurcations of two mode and single mode states, respectively. 
The different types of bifurcations are indicated by SN (saddle 
node), HB (Hopf), PD (period doubling), SNL (saddle node 
of limit cycle), and TR (torus). The striped blue lines TC 
denote transcritical bifurcations between stable single mode 
and stable two mode states. The single mode (stable) lock- 
ing and two mode equilibrium regions are denoted by SL and 
TME, respectively. Only bifurcation lines which affect stable 
states are shown. 



mode power in the region above region (i) in Fig. [TJ As 
we lower the detuning frequency to Alu « 12 GHz, which 
is about twice the relaxation frequency wro , the dynam- 
ical coupling between the uninjected mode and injected 
mode via the carrier density n [Eq. ([3])] becomes relevant, 
and gives rise to a two mode period doubling bifurcation 
(magenta line in Fig. [5]). This bifurcation line forms the 
upper boundary of the dynamical two mode region (i), 
which we experimentally identified in Figs. [T] and [21 

As we further decrease Aw, we cut through a regime 
which is delineated by a torus bifurcation line (red line 
in Fig. [5]) at Aui « 10 GHz. We emphasize that this 
torus bifurcation is a feature the two mode system that 
is not present in the single mode case. Due to this torus 
bifurcation an additional incommensurate low frequency 
of about 300 MHz appears as previously identified in the 
discussion of the characteristic sawtooth antiphase time 
traces of Fig. [T] Note that outside of the boundary of 
the torus bifurcation, the two mode period doubled limit 
cycle is only weakly stable and noise will be able to ex- 
cite low frequency antiphase dynamics. This explains the 
relatively broad region of antiphase dynamics observed 
experimentally in Fig. [2l As we further lower the detun- 
ing frequency to Alu w 7 GHz, the two mode limit cycle 
exchanges stability with a single mode limit cycle in a 
transcritical bifurcation (striped blue hue in Fig. [H]). At 
this point the uninjected mode is completely switched off 
and the limit cycle in the single mode manifold becomes 
stable, forming the lower boundary of the region (i). This 
scenario is in agreement with experimental data of Fig. 



[2l where we observe that the spectral density evolves con- 
tinuously across the lower boundary of region (i). 

As we decrease Alu further, the single mode period 
doubled limit cycle undergoes an inverse period doubling 
bifurcation (dashed magenta line in Fig. [5]) and a sin- 
gle period limit cycle is generated. This is experimen- 
tally verified by the disappearance of the frequency at 
Aw sa 4 GHz, which happens below the lower boundary 
of region (i), i.e. after we have switched to a single mode 
state. 

The generated single period limit cycle, which exper- 
imentally corresponds to the single mode wave mixing 
state, remains stable until we encounter a saddle node 
bifurcation of limit cycles (orange line SNL in Fig. [5]) at 
Aw « 4 GHz. At this point the stable single mode limit 
cycle collides with an unstable one, and both limit cycles 
disappear [2^. The dynamics are then blown out from 
the single mode manifold where the associated bursting 
dynamics explains the brief broadening which is observed 
in the experimental spectrum (Fig. [2]) and in the numer- 
ical spectrum (Fig. [5]). Close to the bifurcation point, 
the remainder of the stable limit cycle still forms a slow 
region on the single mode manifold, and thus the pe- 
riod of the resulting two mode orbit diverges close to the 
bifurcation point. Away from the region of bursting dy- 
namics this new incommensurate frequency results in a 
quasiperiodic orbit on a torus (cf . Fig [8]) . 

Moving away from the SNL bifurcation point to lower 
Aw, the frequency originating at zero increases linearly, 
and gives rise to satellite peaks around the second ma- 
jor frequency close to the relaxation oscillation frequency. 
This explains the striking star-like features in the exper- 
imental (Fig. [2]) and numerical (Fig. [5]) spectra, which 
immediately follow the brief broadening discussed in the 
previous paragraph. Just below the SNL bifurcation, we 
are in a region with two incommensurate frequencies, 
which gives rise to quasiperiodic evolution on a stable 
torus manifold. The torus transforms into a limit cycle 
via an inverse torus bifurcation at Aw w 3 GHz (red line 
in Fig. [9]), and in the power spectrum, the low frequency 
that was introduced by the SNL bifurcation disappears 
again. 

Further lowering Aw below 3 GHz we find that the 
stable two mode limit cycle undergoes a period doubling 
cascade to chaos. In Fig. [HI we only show the first period 
doubling bifurcation (magenta line). At the transition to 
chaotic behaviour the power spectrum broadens dramat- 
ically, which marks the transition from region (ii) to re- 
gion (iii). We leave the chaotic regime at Aw « in a sin- 
gle mode periodic orbit, which undergoes inverse period 
doubling bifurcations (dashed magenta line in Fig. [5]) . 

As in basic single mode laser theory, this limit cycle 
becomes a single mode equilibrium locked state (SL) via 
a Hopf bifurcation (dashed green line in Fig. [5]), which 
disappears in a saddle node bifurcation (dashed black line 
in Fig. [S]). Then the remaining single mode limit cycle 
undergoes again a SNL bifurcation followed by a torus 
bifurcation, as was discussed in the context of region (ii) . 
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This defines the two mode region (iv) of Fig. [TJ 

For even lower Aw we find a two mode equilibrium 
state (TME), which is bounded by a Hopf bifurcation. 
This region corresponds to the experimentally observed 
region below region (iv) in Fig. [1] where the injected 
mode is suppressed. This suppression is a precursor to 
the appearance of a bistability between single-mode in- 
jection locking and a two mode equilibrium at larger neg- 
ative detunings (< —10 GHz) and larger values of the 
injected field strength. We have recently demonstrated 
that this bistability can be the basis of an all-optical 
memory element based on switching of the uninjected 
mode [16i| . Interestingly, there is also a second two mode 
equilibrium region for positive detuning and weak injec- 
tion {0 < Auj <6 GHz, K < 0.005) as shown in Fig. M 
This region will be the subject of more detailed study in 
future work. 



mode laser diode. The device was a specially engineered 
Fabry-Perot laser designed to support two primary modes 
with a large (THz) frequency spacing. At a fixed injected 
field strength, injection in one of the primary modes of 
the device leads to a rich dynamical scenario where dra- 
matic switching between regions of two-mode and single 
mode dynamics was observed as the frequency detuning 
was varied. Using a minimal extension of the single mode 
rate equation model that includes a single equation for 
the carrier density and a general cross and self-saturation 
of the gain, we reproduced experimental results with a 
high degree of accuracy. In particular, we were able to de- 
scribe the dynamical sequence of the total intensity and 
of the individual modes. These dynamics were shown to 
be organised by complex bifurcation structure that re- 
sults from the increased dimensionality of the two-mode 
injected system. 



V. CONCLUSIONS 

We have presented an experimental and theoretical 
study of antiphase dynamics in an optically injected two- 
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